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ABSTRACT 

The Rossby wave instability (RWI) in non-self-gravitating discs can be triggered 
by a bump at a radius ro in the disc surface mass-density (which is proportional to the 
inverse potential vorticity). It gives rise to a growing non-axisymmetric perturbation 
[<x exp(im0), m = 1, 2..] in the vicinity of ro consisting of anticyclonic vortices which 
may facilitate planetesimal growth in protoplanetary discs. Here, we analyze a contin- 
uum of thin disc models ranging from self-gravitating to non-selfgravitating. The key 
quantities determining the stability/instability are: (1) the parameters of the bump 
(or depression) in the disc surface density, (2) the Toomre Q parameter of the disc (a 
non-self-gravitating disc has Q > 1), and (3) the dimensionless azimuthal wavenum- 
ber of the perturbation = mQh/ro, where h is the half-thickness of the disc. For 
discs stable to axisymmetric perturbations (Q > 1), the self-gravity has a significant 
role for k^ < ir/2 or m < (ir /2)(ro / tyQ^ 1 ; instability may occur for a depression or 
groove in the surface density if Q < 2. For k^ > tt/2 the self-gravity is not important, 
and instability may occur at a bump in the surface density. Thus, for all mode num- 
bers m ^ 1, the self-gravity is unimportant for Q > (%/2)(ro/h). We suggest that the 
self-gravity be included in simulations for cases where Q < (ro/h). 
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1 INTRODUCTION 

The theory of the Rossby wave instability (RWI) was de- 
veloped by Lovelace et al. (1999) and Li et al. (2000) for 
thin accretion discs with negligible self-gravity and ear- 
lier by Lovelace & Hohlfeld (1978) for thin disc galaxies 
where the self-gravity may or may not be important. In 
the first case the instability can occur if there is a bump 
as a function of radius in the inverse potential vorticity 
F(r) oc Efi[(V x u) ■ z]" 1 at some radius ro, where £ is 
the surface mass density, the angular velocity, and u the 
flow velocity of the disc. Such a bump could arise at the ra- 
dial boundary of the "dead zone" (Varniere & Tagger 2006; 
Lyra et al. 2009; Crespe et al. 2011 ). This zone can arise 
from the suppression of the magnetorotational instability by 
low ionization inside the disc (Gammie 1996). In turn, the 
bump in F(r) can give rise to the exponential growth of non- 
axisymmetric perturbations in the vicinity of ro consisting 
of anticyclonic vortices. Such vortices can act to concen- 
trate dust grains in their centres and thereby accelerate the 
formation of macroscopic planetesimals (Barge & Somme- 
ria 1995; Tanga et al. 1996; Bracco et al. 1999; Godon & 



Livio 2000; Johansen et al. 2004; Heng & Kenyon 2010). 
Two-dimensional hydrodynamic simulations of the RWI in- 
stability in discs were done by Li et al. (2001) and recently in 
three dimensions (e.g., Meheut et al. 2010, 2012a). The the- 
ory of the instability in three dimensions has been studied 
by Meheut et al. (2012b) and Lin (2012). The instability is 
also predicted to occur in strongly non-Keplerian discs with 
regions where dQ/dr > which exist around rotating mag- 
netized stars (Lovelace et al. 2009). In the RWI the wave is 
radially trapped within the disc which is different from the 
Papaloizou and Pringle (1984; 1985) instability where the 
wave is trapped between the inner and outer radii of a disc 
or torus. The RWI has an important role in the accretion- 
ejection instability of discs proposed by Tagger and collabo- 
rators (e.g., Tagger & Varniere 2006; Tagger & Melia 2006; 
Tagger & Pellat 1999). 

The RWI in disc galaxies consisting mainly of stars has 
been studied using JV— body simulations (Sellwood & Kahn 
1991; and Sellwood 2012). When the self-gravity of the disc 
is important, instability occurs at radii where F(r) has a 
minimum, that is, a groove (Lovelace & Hohlfeld 1978; Sell- 
wood &Kahn 1991). We do not consider the instability of 
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self-gravitating gas discs for conditions where the gas cool- 
ing is important (Gammie 2001). 

Saturation of the exponential growth of the instability 
is predicted to occur when the time-scale for a fluid par- 
ticle to orbit the center of a vortex is comparable to the 
wave's growth time-scale (Lovelace et al. 2009). Simulations 
by Meheut et al. 2012c) support this conclusion. 

Section 2 develops the theory, applies it to radially local- 
ized perturbations, and considers axisymmetric and nonax- 
isymmetric modes. Section 3 gives sample results for cases 
where self-gravity is important and where it is negligible, 
and Sec. 4 gives a numerical example relevant to forming 
protostar. Conclusions are given in Sec. 4. 



2 THEORY 

We consider the stability of of a thin self-gravitating disc 
of equilibrium surface mass density E(r) with a (r, cf>, z) 
coordinate system. The equilibrium has the flow velocity 
u = M^(r)0 = rQ(r)(p . That is, the accretion velocity 
u r and the vertical velocity u z are assumed negligible com- 
pared with v,(j>. The equilibrium flow satisfies — Erf2 2 = 
—dP/dr — EV3>, where P the vertically integrated pressure 
and $ the gravitational potential. This potential is give by 
V 2 $ = 4tyGT,8(z), where G is the gravitational constant. 

The perturbed quantities are: the density, E = E + 
<5E(r, 4>, t); the pressure is P = P + 8P(r, <f>, t); the flow ve- 
locity is u = u + 5u(r, 4>, t) with 8u = (5u r , 8u$, 0). The 
equations for the perturbed flow are 



From equation (lc) and (Id), we have 



PS 
~Dt 



+ E V • u = , 



(la) 



SP- 



c 2 c5E- 



AuLs 



5u r 



Here, c a = (dP/dT,)^ 2 is the effective sound speed in the 
disc and Lg 1 = ■y~ 1 dln(S)/dr with L.s the length-scale of 
the entropy S = P/E 7 variation in the disc. To simplify the 
subsequent calculations we consider the homentropic case 
where Ls — > oo. 

The perturbation of the gravitational potential is give 

by 



V 2 <5$ = 4irG5T,5(z) . 
Equations (3) can be solved to give 
~Alj 95* 



Eou r — iF 



ESuj, = F 



fi dr 
Q. 2 dSt> 



2D. 2 dr 



2k^,5<5> 



k,— 5* 



where 



_ 2 <5E 



is an effective potential and 



F = 



En 



Q 2 - (Au) 2 



(4) 

(5a) 
(56) 

(5c) 

(5d) 



is the function identified by Lovelace and Hohlfeld (1978). 
Substituting equations (5) into (2) we obtain 
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where D/Dt = d/dt + u • V, and where S = P/(E) 7 is the 
entropy of the disc matter. 

We consider perturbations ~ f(r)exp(im(j> — iijjt), where 



2.1 Radially Localized Modes 

Here we consider radially localized modes in the sense that 
perturbation extends over a radial region Ar with (Ar) 2 <g 
r 2 centred at r . With this in mind we rewrite equation (6) 

as 

2 2 ST, d 2 0* 2 

(fir "AW )— = -p--fc^tf 



gular frequency. For free perturbations cj = uj r + iuii and for 
the growing modes of interest > 0. From equation (la), 
we have 



where 



iAuj <5E = V • (E <5u) , 



Au(r) ee to — mfi(r) 



(2) 



and — u^/r. 

From equation (lb) we have 

i Au8u r + 2Q,8u 4 , = i ^ -^T^r- + Vo$ , (3a) 
T, dr T, 2 dr 

q 2 §p 
iAcoSu,/, — ^rfrSur = ik<f,— + ik<f,S& (36) 

Here, f2 r = [r~ 3 d(r 4 f2 2 )/dr]3 is the radial epicyclic fre- 
quency, and k$ = m/r is the azimuthal wavenumber. For 
an approximately Keplerian disc, S~l r « ^. 



d In F dS9 2k lj> £ld\nF ^ 
dr dr Au) dr 



where F = rF/Q. 
We let 



/dr 
— (<5E, «5*)expHMr-r„)] , 

where k r is the radial wavenumber of the perturbation. The 
radial Fourier transform of equation (7) gives 



+ 



lnFd(5* 2k 4 ,Q,d\nF 



dr dr Auj dr 



5* 



(8) 



Here, we have neglected the radial variation of E and Q r , and 
that of Au) 2 in comparison with fl 2 . Also, k = k r i + k^,(p 
and (..)fc r denotes the Fourier transform. 

Assuming \k r \ ~ (Ar) -1 3> ro, the WKBJ solution of 
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Figure 1. The horizontal lines mark out the regions of 'bump' 
and 'groove' instability in a self-gravitating disc for m = 3 and 
the ring instability for m = 0. Here, A c = 2ir/k c and A^ = 2-rr/k^. 
For the Q > 1 part of the figure we have used equations (12) and 
(15) with AF = 0.2 for the bump and AF = -0.2 for the groove, 
Ax = 0.04, and q = 3/2. 



equation (4) gives 8&k r = — 27rG5£fc r /|k|. Therefore, from 
equation (5c) we obtain 



where 



<5* fc „ = (l--^U — 



ttGE 



(9a) 



(96) 



is a characteristic wavenumber. 



2.1.1 Axisymmetric Modes 



For axisymmetric perturbations (k^ = = m) of smooth 
disc (where the term oc dhiF/dr can be neglected in equa- 
tion (8) ), equations (8) and (9a) give 



(10) 



(11) 



(Safranov 1960; Toomre 1964). The minimum of [ui(k r )] 2 oc- 
curs at k r — k c where uj 2 = Q, 2 — k 2 c 2 . Therefore, if we define 
fc c * = fl r /cs, then for k c < k c * the minimum of a; 2 is positive 
and the axisymmetric perturbations are stable. Conversely 
for k c > fc c * the perturbations are unstable. With 

k c 7rG£ 

the axisymmetric perturbations are stable (unstable) for 
Q > 1 (Q < 1) (Toomre 1964). The minimum value of the 
squared frequency is uj(k r ) 2 = (fc c ,c s ) 2 (l — Q~ 2 ). Figure 1 
shows the unstable region for a ring, Q < 1. 

2.1.2 N on- Axisymmetric Modes 

For the non-axisymmetric perturbations (m = 1,2, .., — 
m/r) we assume that (Aw) 2 <C fi 2 in equation (8). Multi- 
plying equation (8) by (1 — 2fc c /|k|) and doing the inverse 
Fourier transform gives 

Q. dlnF , , . 

— 8^(x ) . (12a) 



2 f'°° 

5V (x) = - / dx'K(x ■ 



Au(x') dx 



This represents an integral equation for S^(x). The term 
in equation (8) proportional to d In F/dr has been dropped. 
It is negligible compared with the term involving d In F/dr 
because |mfl/ ' Au>\ 3> 1 for the considered modes. The kernel 



K(x) = - 

7T 



cos(mxy) 



dy 



- 1 



1 + V 2 + $ 



(126) 



Here, x = (r — ro)/ro with i 2 < 1, !1 = Q(ro), Q is given in 
equation (11), and k^ = k^/k c . 

The half-thickness of the disc h is c s /£l r so that 



mQ 



(13) 



where we assume (h/r) 2 <C 1. Note that A c = 2ix/k c = 
2nhQ. The Toomre (1964) characteristic wavelength is de- 
fined as A* = 47rGE/n 2 = Ah/Q. Thus A c /A* = ttQ 2 /2. 
In parallel with equations (12) note that 

2 *s(s) 2 /■« , , n 

c, — — — — — I dx G(x — x ) — — t- , 

s E 7T 7.^ k 1 Aw{x') dx 

(14a) 

Here, 



dlnF rT , , 
— ) 



G(x) 



cos(ma;3/) 



l + y 2 + 



<7- 



2\/l+y' 



(146) 



is a second kernel. 

We solve equation (12) for the case where ln[F(a;)] has 
a square bump or groove. That is, 



ln[F(x)] = AFH{x + Ax/2)H(Ax/2 - x) 



(15) 



where H is the Heaviside step function (+1 for positive ar- 
gument and zero for a negative argument), and Ax is the 
fractional radial width of the bump (AF > 0) or groove 
(AF < 0). This distribution of InF is known to be useful 
for analytic calculations of the non-axisymmetric instability 
of discs (Sellwood & Kahn 1991; Umurhan 2010). 
Substituting equation (15) into (12) gives 



2 „ QAF . 

7T Auj- 



2 „ QAF . T 
■k flu_ 



2 O.AF . _ 

7T ALU+ 



2 „ Q.AF . 
-K — o*+ 



(16a) 



(166) 



where Ko = F'(O), Kax = K(Ax) and Aui± — uj— mf2[ro(l± 
Ax/2)] ~ uj—mQiqm AxQ/2 for the g = — din fi/dlnr with 
q — 3/2 for a Keplerian disc. Equations (16) give a quadratic 
in uj which for the assumed symmetry of the bump or groove 
about x — gives either (w — mQ) 2 > - stable motion - or 
(uj — mQ) 2 < - unstable motion. 

Evidently an arbitrary ln[F(r)] can be approximated by 
a staircase-like function with the result that equation (12) 
leads to an n x n matrix equation in place of equation (16) 
with n the number of steps. 



3 SAMPLE RESULTS 

The main parameters of importance for instability are: (1) 
The parameters of the bump or groove (AF and Ax). (2) 
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Figure 2. Growth rates of axisymmetric perturbations (m = 0) 
for Q < 1 and non-axisymmetric perturbations (m = 3, 5) for 
Q > 1 for a groove (AF < 0) in the F(r) profile, k^ = A c /A^ = 
0.5, and q = 3/2. Here, AF C is the threshold value of the depth of 
the groove AF needed for instability. Also, Ti is proportional to 
the radial angular momentum flux across the centre of the groove 
as discussed in the text. For the Q > 1 part of this figure we have 
used equations (12) and (15) with AF = —0.2 for the groove and 
Ax = 0.04. 



The parameters of the disc Q and of the perturbation m 
and k<f, — k^/kc. Figure 1 gives a qualitative picture of the 
regions of bump, groove, and ring instability in a disc as a 
function of Q and k^, = k^fk c . The needed kernel in equation 
(12b) is accurately evaluated numerically using a 4000 point 
integration. (Our attempt to analytically integrate (12b) was 
thwarted by the required branch cuts at the points y = ±i.) 
The transition from groove to bump instability occurs when 
K(0) goes through zero as a function of k<j>, with K(0) > 
for small k<j>. We find that this occurs at » 1.55. This 
agrees with the conclusion of Lovelace and Hohlfeld (1978) 
that the transition is at kj, = n/2 w 1.57. 

Figure 2 shows the growth rates as a function of Q for 
the ring mode (m = and Q < 1) and the groove mode 
(m = 3, 5, Q > 1, and k$ = 0.5). Also shown in the fig- 
ure is AF C which is the minimum groove depth needed for 
instability. Further, Tt is proportional to the radial angu- 
lar momentum flux across the x = surface. Specifically, 
Tt oc r J d(f>?R.(r 2 QST,* 8u r + T,r8u^8u r ), where 5ft denotes 
the real part. 

Generally, for both unstable grooves and bumps the 
growth rate increases as \AF\ increases above a threshold 
value. Also, the growth rate decreases as the width of the 
groove or bump Aa; increases. This supports the interpre- 
tation of Umurhan (2010) that the instability of a bump 
arises from interaction "edge waves" on the surfaces where 
F changes rapidly. 

Unstable grooves: Figure 3 shows the density pertur- 
bation <5£(r, cj>) and the arrows represent the flow perturba- 
tion (Su r , 5u$) for an unstable groove with AF = —0.2, 
and Q = 1.3, m = 3, and k<f, = 0.5. Notice that the an- 
ticyclonic motion of the vortex about the minimum of the 
density perturbation. 

Unstable bumps: Figure 4 shows the density pertur- 
bation <5£(r, <j>) and the arrows represent the flow perturba- 
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Figure 3. The background gray-scale (white=low, dark=high) 
and the contours represent the density perturbation <5£(r, <f>), and 
the arrows represent the flow perturbation (8u r , Su^) for the 
case of an unstable groove, AF = —0.2, and Q = 1.3, m = 3, 
= 0.5, and q = 3/2. The growth rate is £jj/f2 = 0.152. The 
larger arrows near the top and bottom of the figure have been 
put in to make clear the direction of the circulation which is 
anticyclonic. In this case this circulation is around a minimum of 
the density perturbation. 
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Figure 4. The background gray-scale (white=low, dark=high) 
and the contours represent the density perturbation <5S(r, <j>), and 
the arrows represent the flow perturbation (5u r , Su^) for the case 
of an unstable bump, AF = 0.2, and Q = 1.3, m = 3, k^ = 2, and 
q = 3/2. The growth rate is ujj/fi = 0.0844. The larger arrows 
near the top and bottom of the figure have been put in to make 
clear the direction of the circulation which is anticyclonic. In this 
case this circulation is approximately around a maximum of the 
density perturbation. 



tion (8u r , 5u<f,) for an unstable bump with AF = 0.2, and 
Q — 1.3, m = 3, and k$ = 2. Notice that the center of the 
anticyclonic motion of the vortex is offset by a small angle 
from the maximum of the density perturbation. The angular 
momentum flux across x — is found to be always outward 
for unstable bumps. 

Figure 5 shows a further case of an unstable bump with 
AF = 0.2, and Q — 2, m — 3, and fc^ = 2. The center of the 
anticyclonic motion of the vortex is again offset by a small 
angle from the maximum of the density perturbation. 



4 NUMERICAL EXAMPLES 

We first consider an accretion disc around a solar mass star 
with accretion rate M = 10 -8 Moyr _1 . For an a— viscosity 
disc with a = 5 x 10 -3 , the surface mass density is E = 
89.3 gem -2 at r = 10AU assuming the half-thickness of the 
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Figure 5. The background gray-scale (white=low, dark=high) 
and the contours represent the density perturbation <5E(r, <f>), and 
the arrows represent the flow perturbation (Su r , Su^) for the case 
of an unstable bump, AF = 0.2, and Q = 2, m = 3, = 2, and 
q = 3/2. The growth rate is u>i/Q = 0.0661. The larger arrows 
near the top and bottom of the figure have been put in to make 
clear the direction of the circulation which is anticyclonic. In this 
case this circulation is approximately around a maximum of the 
density perturbation. 



clonic vortices have density enhancements. The growth rates 
for both the groove and bump instability are < 10% of the 
orbital angular frequency for bumps fractional amplitude 
<20%. 

For fixed radius ro and sound speed c s , the regime of 
instability - groove or bump - is set by the Q value which 
is proportional to M/a with a the dimensionless viscos- 
ity. For Q < 2 and m < (tv /2)(ro/h)Q~ 1 the groove in- 
stability may occur. On the other hand for Q ^> 1 and 
m > (n/2)(ro/h)Q~ 1 the bump instability may occur. The 
estimates of Sec. 4 indicate that the two regimes are not 
widely separated in terms of M and a. For this reason wc 
suggest that simulations of Rossby vortices include the self- 
gravity for cases where Q < r /h. 



ACKNOWLEDGEMENTS 

We thank the referee for helpful criticism. RVEL was sup- 
ported in part by NASA grant NNX11AF33G. 



disc is h — 0.04r or c s /vk = 0.04, where c s is the isothermal 
sound speed in the disc and Vk = 9.4 kms -1 is the Keplerian 
velocity. We have used M — 27rrE« r with u r — ac^/vx- 
For comparison, the surface mass-density of a Jupiter mass 
spread over a radial extent 2h at 10AU is about 170 gcm~ . 
From equations (9b) and (10), k c = 1.97/r and Q = 12.7. 
For fixed c s /vk, notice that Q oc M/a. In equation (13) wc 
have kj, — 0.04m<2 = 0.508m. For k<f, < it/2 or m < 3 there 
is no instability because Q> 1. We have a bump instability 
for m > 3. Clearly, all of the modes m ^ 1 are in the bump 
regime if Q > ■nr/(2h). The mass of the disc is relatively 
small, ~ 7rr 2 E/M Q = {c s /v K )/Q = 0.0031, so that it is 
approximately Keplerian. 

As a second case, consider a larger accretion rate M = 
3 x lO _8 M0yr _1 and a smaller viscosity a = 0.002 but with 
the other quantities (r, c s ) the same as above. Then we 
find that Q = 1.69 and fc = 0.0676m. For fc < tt/2 or 
rn < 23 there can be a groove instability because Q is not 
much larger than unity (see Fig. 1). At the same time we 
expect the bump instability for m > 23 to be unimportant 
because the azimuthal wavelengths are comparable to the 
disc half-thickness. 



5 CONCLUSIONS 

We have analyzed the Rossby wave instability in a con- 
tinuum of thin disc models ranging from self-gravitating 
to non-selfgravitating. The important quantities determin- 
ing the stability/instability are: (1) the parameters of the 
bump (or depression) in the inverse potential vorticity F(r) 
at ro, (2) the Toomre Q parameter of the disc, and (3) 
the dimensionless azimuthal wavenumber of the perturba- 
tion k^ = mQh/ro, where h is the half-thickness of the 
disc and m = 1, 2, .. is the azimuthal mode number. For 
k^ < tt/2 and Q < 2 instability may occur for the case of 
a groove in the surface-density. In this case the centres of 
the anticyclonic vortices have reduced density. For fc</> < 7t/2 
Q > 2, both bumps and grooves are stable. For k^ > tv/2, 
the bumps may be unstable and the centres of the anticy- 
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